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1.0 INTRODUCTION 


In aircraft turbine engine hot section components, cyclic thermal 
stresses are the most important damage mechanisms. Consequently, an accurate 
and reliable prediction of thermal loads is essential to improved durability. 

To achieve this goal, considerable effort has been expended over the past 20 
years to (1) acquire engine temperature test data, and (2) develop accurate, 
reliable, efficient computer codes to predict steady-state and transient temp- 
eratures and for calculating elastic and inelastic cyclic stresses and strains 
in hot section components. There is a need for continued development of these 
codes because the availability of more accurate analysis techniques for complex 
configurations has enabled engine designers to use more sophisticated designs 
to achieve higher cycle efficiency and reduce weight. 

It became apparent m recent years that a serious problem existed when 
interfacing output temperatures and temperature gradients from the heat 
transfer codes with the input to the stress analysis codes. Due to the growth 
m computer capacity and speed, and with the development of input preproces- 
sors and output postprocessors, the analysis of components using hundreds or 
thousands of nodes in the heat transfer and stress models became economical 
and routine. This has aggravated the problem of manually transferring output 
temperatures from heat transfer nodes to stress analysis input to the point 
where the engineering effort required is comparable to that required for the 
remainder of the analysis. Furthermore, a considerable amount of approxima- 
tion has been introduced in an effort to accelerate the process. This tends 
to introduce errors into the temperature data, thereby negating the improved 
accuracy in the temperature distribution achieved through use of a finer 
mesh. 


The objective of the Burner Liner/Thermal Structural Load Modeling pro- 
gram was to develop a transfer module that would accurately and efficiently 
transfer three-dimensional temperature information from heat transfer analysis 
codes to stress analysis codes. This objective has been satisfied by the 
development of the computer code TRANCITS (TRansfer ANalysis Code to ^Inter- 
face Thermal/Structural Problems). 

This report is a preliminary user's manual for the TRANCITS program. 
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2.0 GENERAL OVERVIEW 


TRANCITS is a computer code developed by the Aircraft Engine Business 
Group of the General Electric Co. in conjunction with the NASA-Lewis Research 
Center to address the problem of transferring three-dimensional temperature 
information from numerical heat transfer analyses to numerical stress 
analyses. The advanced features of TRANCITS include: 

• The ability to transfer temperatures between dissimilar mesh 
densities 

• The ability to handle finite difference and finite element heat 
transfer codes 

• Extensive user controls on the interfacing procedure 

• Automatic interface to SINDA, MARC, and NASTRAN 

• Internal checks on the accuracy of the transfer 

• The ability to handle stress points slightly outside the heat 
transfer model 

• Accurate and efficient transfer of temperatures 

• Modular program architecture to facilitate maintenance and 
enhancements . 

Figure 1 is a schematic of how TRANCITS is used in the thermal interfacing 
procedure . 

It should be noted that although the interfaces between SINDA, MARC, and 
NASTRAN codes are Hard Wired into the TRANCITS system, the Neutral Input and 
Neutral Output files allow TRANCITS to be used with any combination of heat 
transfer codes and stress analysis codes. 
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3.0 CURRENT LIMITATIONS 


The following are the current limitations for the first release of 
TRANCITS (Version 1.0). 

1. Heat Transfer Model Limitations 

* 

• Maximum number of heat transfer limitations = 5000 

• Maximum heat transfer element name = 5000 

• Maximum number of transient solutions = 500 

• Maximum number of heat transfer nodes = 5000 

• Maximum heat transfer node name = 5000 

• Heat transfer element library - 8-noded 
solid element 


2. Stress Analysis Model 

• Maximum number of stress nodes 

• Maximum stress node name 

• Maximum number of stress elements 

• Maximum stress element name 


= 5000 

= 5000 

= Unlimited 
= Unlimited 


Stress element library - any three-dimensional 
solid element 


All of the numerical limits can be easily increased by increasing the amount 
of temporary storage allocated to the program. 
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4.0 DESCRIPTION OF INPUT 


TRANCITS can accept many different input forms, including the output of 
the SINDA and MARC programs directly. The choice of the input form depends 
on which heat transfer code was used to produce the temperatures. For heat 
transfer codes other than MARC and SINDA, a neutral heat transfer input file 
can be used to input the required data into TRANCITS. 

Both finite difference and finite element codes can be used, but in 
either case the geometry of the heat transfer mode and the computed 
temperatures from the heat transfer analysis must be input. In addition to 
this information, a user control file must be entered. This file contains 
data that controls the interface procedure. All of these input files are 
described in the following sections. 


4.1 NEUTRAL HEAT TRANSFER INPUT FILE 


This file is either input or created internally. It contains all the 
geometry and temperature information for the heat transfer model. All 
processing of this information in TRANCITS is done using this file. If the 
heat transfer code used was SINDA or MARC, the output of these two programs 
is input into TRANCITS and the Neutral Heat Transfer Input file is created 
internally. On the other hand, if the heat transfer code used was some other 
program, the output of this program could be externally converted to the 
Neutral Heat Transfer Input form and then interfaced by using TRANCITS. The 
use of the Neutral Heat Transfer Input file allows TRANCITS to be used with 
any heat transfer analyzer. 

This file is composed of six sets, or partitions, of data. The six 

are : 

TITLE 
DIRECTORY 
NODAL 

ELEMENT CONNECTIVITY 
ELEMENT FACE 
TEMPERATURE 

Title Partition 

The title partition allows the input file to have up to four 80- 
character titles. These titles appear on the first four records of the 
input file m the form: 

TITLE 1 RECORD 1 

TITLE 2 RECORD 2 

TITLE 3 RECORD 3 

TITLE 4 RECORD 4 


PARTITION 

PARTITION 

PARTITION 

PARTITION 

PARTITION 

PARTITION 


5 



Directory Partition 


The directory partition contains count and control information about the 
input. It is of the form: 


NSHIFT, NUMN, NUME, NUMT, NTRECS , ECODE, ACODE 


where 


NSHIFT 

NUMN 

NUME 


No. of records in the directory partition - 1 
(currently NSHIFT = 0) 

No. of nodes in the heat transfer model 
No. elements in the heat transfer model 


NUMT 

NTRECS 


ECODE 


No. of transient solutions on the input file 

No. of heat transfer elements with temperature information 
(usually NTRECS = NUME for ACODE = 1 or NTRECS = NUMN for 
ACODE = 2) 


Heat transfer element type code 


1- 4 noded 

2- 8 noded 

11- 8 noded 

12- 16 noded 

13- 20 noded 


2D element 
2D element 
3D element 
3D element 
3D element 


ACODE = Heat transfer analysis type code 

1 - Finite difference 

2 - Finite element 


Nodal Partition 


The next NUMN records compose the nodal directory. It contains the 
rectangular cartesian coordinates for every node in the heat transfer model. 
The form of this directory is: 

Node X Y Z 

. . . . NUMN Records 

where NODE = Node Name 
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X 

Y = Coordinates 
Z 


Element Connectivity Partition 

The next NUME records in the file compose the element connectivity parti- 
tion. It contains all the heat transfer element connectivity information in 
the form: 

NE IE(L) , L = 1, NUNPE 
. . NUME Records 


where 

NE = Element name 
IE = Connectivity array 
NUNPE = No. of nodes each element 


Element Face Partition 


This partition is omitted if the heat transfer analysis type was finite 
element (ACODE = 2) . If ACODE = 1 (finite difference) , this partition 
contains NUME records of the form: 

NE IFACP 


NUME Records 


where 

NE = Element name 
IFACP = Packed element face number 

IFACP is a six-digit integer that relates the finite difference face labels 
to the sides of the heat transfer element. It is required when the finite 
difference element has the ability to produce face center temperatures as 
well as element centroidal temperatures. This input allows TRANCITS to use 
all temperature information produced by the heat transfer element in the 
interfacing procedure. Figure 2 shows how IFACP is defined for an eight- 
noded solid heat transfer element. 

SPECIAL NOTE : When the heat transfer element being used produces only 

centroidal temperatures, IFACP should be set to 123456. Therefore, for SINDA 
applications IFACP = 123456. 
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Temperature Partition 


This partition contains all temperature information from the heat 
transfer analysis. It can appear in two forms, depending on the type of heat 
transfer analysis. 


- FINITE DIFFERENCE TEMPERATURE PARTITION - 
(ACODE = 1) 

The form is: 

TIMEN 

Solution NE TMP(L) , L = 1, NTPE 

1 

(NTRES + 1) Records 


Solution 

2 


TIMEN 

NE TMP(L) , L = 1, NTPE 

. . (NTRECS + 1) Records 


Solution 

NUMT 


(NTRECS + 1) Records 


- FINITE ELEMENT TEMPERATURE DIRECTORY - 
(ACODE = 2) 


The form is: 


Solution 

1 


TIMEN 

NODE TP 


(NTRECS + 1) Records 


Solution 

2 


TIMEN 

NODE TP 


(NTRECS + 1) Records 
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Solution 

NUMT 


(NTRECS + 1) Records 


where 


Time N 


Time value of transient solution 

NE 

- 

Element name 

TMP 

= 

Temperature array 

NTPE 

= 

No. of temperatures each elemet 


= 

5 for 2D 



7 for 3D 

NODE 


Node name 

TP 

= 

Nodal temperature 


An example of this file is given m Appendix B. 


4.2 SINDA TEMPERATURE INPUT FILE 


This file is necessary if the heat transfer analysis code was SINDA. 
The users of the SINDA program have the capability of formatting the out- 
put in many ways. This formatting is accomplished by including Fortran 
subroutines in the SINDA job stream to control the output. An example of 
these output routines for SINDA is given in Appendix C. The SINDA output 
produced by these routines is defined as the SINDA Temperature Input File. 
The form of this file is: 


NUME 


NE1 NE2 .... NE10 

NE11 NE12 .... NE20 

NUME 

Elements 


TIME1 

T1 T2 

T7 T8 

TIME2 

T1 T2 

T7 T8 


T6 

T12 

NUME 

Temperatures 


T6 

T12 

NUME 

Temperatures 


etc . 
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where 


NUME = No. of heat transfer elements 

NXX = Element name 

TXX = Temperature for the element (in the same order as elements) 
TIMEXX = Time value for this transient solution 

This input form is the required input into TRANCITS. It can be obtained 
through SINDA by using the routines given m Appendix C or by external 
formatting of other SINDA output files. An example of the SINDA temperature 
input file is given in Appendix B. 


4.3 SINDA GEOMETRY INPUT 

This input is needed only if the heat transfer analysis code was SINDA. 
The SINDA output files do not contain geometry information; therefore, the 
SINDA GEOMETRY FILE must be input separately into TRANCITS. This file 
contains the coordinate information for all the nodes in the heat transfer 
mode as well as the connectivity of all the heat transfer elements. The form 
of the file is: 

TITLE 

NUMNOD NUMELM NUNPE 

NODE X Y Z 

. . . . NUMNOD Records 

NE, (IE(L), L = 1, NUNPE), IFACP 

. . . NUMELM Records 


where: 

TITLE = 80-Character title 

NUMNOD = Number of nodes m the heat transfer model 
NUME = Number of elements in the heat transfer model 

NUNPE = Number of nodes each element 

NODE = Node name 

X,Y,Z = Rectangular cartesian coordinates of heat transfer nodes 
NE = Element name 

IE = Connectivity array of heat transfer element 

IFACP = 123456 for SINDA 

All data in this file are free format, but each node and element record must 
appear on a separate line of the file. An example of this file is given 
in Appendix B. 
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4.4 MARC HEAT TRANSFER INPUT 


Heat transfer information from a MARC posttape can be input to the 
TRANCITS program. The MARC posttape can be in either formatted or sequential 
binary form; however, the following limitations apply: 

♦ Coordinate and connectivity data must be included on the posttape. 
This option is turned on with the postmodel definition card. 

♦ Only eight-noded heat transfer solids (Element Type 23) are read 
from the tape and used in TRANCITS. 

♦ Rezoning (that is, changes m coordinate and connectivity data) is 
not permitted. 

The contents of the posttape are assumed to be as described in Appendix 
A. 


MARC heat transfer input is selected by setting the variable IHTIC in 
the namelist input file to 2 for a formatted posttape or 3 for a sequential 
binary posttape. 

An example of the formatted MARC postfile is given in Appendix B. 


4.5 STRESS MODEL COORDINATE INPUT FILE 


This file is not optional. It must be input. It contains the 
coordinate information for each stress node in the stress analysis model. 
The form of the file is: 

Node X Y Z 


where: Node is the STRESS NODE LABEL 

X Rectangular 

Y Cartesian coordinates of the stress node 

Z 

For an example see Appendix B. 

The data are in free formatted form, but each nodal record must appear 
on a separate line of the file. 

Every stress node in this file will be interfaced with the appropriate 
heat transfer element m the heat transfer model. 
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4.6 STRESS MODEL ELEMENT INPUT 


This file is optional. It must be input only if MARC Gauss point 
temperatures or MARC centroid temperatures are requested. It contains the 
connectivity information for the stress model elements. The form of the file 
is : 

NUMETP 

NETYP NUME 

Set 1 NE (IE(L) , L = 1, NUNPE) 

No. of records = No. of elements +1 
. . of this type 

NETYP NUME 

Set 2 NE (IE(L), L = 1, NUNPE) 

. . No. of records = No. of elements +1 

. . of this type 

Set 3 


Set 

NUMETP 

where 

NUMETP = No. of sets of element types 
NETYP = Element type code 
Code Element 


1 

4 

noded 

2D 

element 

2 

8 

noded 

2D 

element 

11 

8 

noded 

3D 

element 

12 

16 

noded 

3D 

element 

13 

20 

noded 

3D 

element 


NE = Element name 

IE(L) = Element connectivity 

NUNPE = No. of nodes each element 

These data are free formatted, but each element record must appear on a 
separate line of the file. See Appendix B for an example of this file type. 
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If this file is not input and if MARC temperatures are requested, only 
the temperatures of the stress points will be output. 


4.7 USER CONTROL FILE OUTPUT 


This file is not optional. It must be input. It contains input varia- 
bles that control the interfacing procedure. This file is m namelist format 
which allows the user to input only the variables related to the problem of 
interest. All other variables will default to internal values set by 
TRANCITS. The form of the file is: 

(2) (starts in Column 2) 

$ZMAF 

Variables 

$ 


All variables currently available and their defaults are described m the 
following tabulation. 


Variable Valid Default 

Name ^Values Value 


Description 


I TRAN 0 

1 


0 Do not transform heat transfer coordinates. 

Transform heat transfer coordinates. 


XOFF 

YOFF 

ZOFF 


0 If ITRAN =1, add the values of XOFF, YOFF, and 

0 ZOFF to the coordinates of the heat transfer 

0 nodes , 


ISL 0 

1 


XMIN 

YMIN 

ZMIN 

XMAX 

YMAX 

ZMAX 


NUMZAP/ 

NZAP 


0 Do not section the heat transfer model. 

Section the heat transfer model. 

-10 5 If ISL = 1, these minimum and maximum 

-10 s coordinates determine which heat transfer 
-10 s elements will be considered for the inter- 

-10 5 facing procedure. Note, if any node of a 

-10 5 heat transfer element is outside the window 

-10 s defined by these mm's and max's, the heat 

transfer element is not used. These coordi- 
nates are m the transferred heat transfer 
coordinate system. 

0 This input is used to eliminate unwanted heat 

transfer elements. A shorthand input is avail- 
able (such as the input NUMZAP/NZAP = 1, -7, 

-2, 12, 101, -104, 203, -218, -3) which would 
eliminate heat transfer elements 1, 3, 5, 7, 

12, 101, 102, 103, 104, 203, 206, 209, 212, 215, 
and 218. 
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Variable Valid Default 

Name Values Value 


Description 


NTSTP/ 

TIMD 


NTMDO 


NSKIP 


TIMESP 


IHTIC 


ISPOC 


0 

1 

2 

3 


0 

1 

2 

3 


Note, a maximum of 200 numbers may be entered 
using the NZAP array. 

-1 This input is used to select the desired tran- 
sient solutions from the heat transfer input 
file. If it is omitted or entered as NTSTP/ 

TIMD = -1, all solutions will be interfaced 
(such as the input NTSTP/TIMD = 10, 20, 31, 
would select three solutions 10 s, 20 s, 31 s). 

A maximum of 500 solutions may be selected. 

999 This input allows the first NTMDO solutions to 
be selected. For example, NTMDO = 5 would 
select the first 5 solutions. 

0 This input allows the first NSKIP solutions to 

be skipped or not processed. For example, NSKIP 
= 50, NTMDO = 3, would skip the first 50 solu- 
tions and process the next 3. 

IE-5 This input is the tolerance between the desired 
time value and the time value on the heat trans- 
fer input. If TIMESP = 0.01 and TIMD = 10, any 
time values between 9.99 and 10.01 would be 
selected. 

Heat Transfer Input Code 

0 Input is Neutral Heat Transfer Input file 

Input is SINDA temperature file 
Input is MARC temperature file (formatted) 

Input is MARC temperature file 
(binary sequential) 

Stress Program Output Code 

0 Output is Neutral Temperature file 

Output is formatted for NASTRAN 
Output is formatted for MARC 
(element controid) 

Output is formatted for MARC (element Gauss 
point) 

2X2X2 


Note that any or all of these varaiables can be input using the User 
Control Input file. At least one variable must be input. It is recom- 
mended that if all defaults are desired, the input file should be: 
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$ZMAF 

NTMDO = 999 , 

$ 

An example of the file is given in Appendix B. 
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5.0 OUTPUT DESCRIPTION 


TRANCITS produces a variety of different output files depending on which 
stress analysis output code has been selected. Two files are always produced, 
the hard-copy file and the neutral temperature output file. The two optional 
output files are the NASTRAN temperature file and the MARC temperature file. 
The following sections describe all of these various outputs. 


5.1 NEUTRAL TEMPERATURE OUTPUT FILE 

The Neutral Temperature Output file is the primary output of TRANCITS. 

It contains the interfaced temperature for all of the stress points processed 
for every transient solution selected. The form of the file is: 

Formats 

2X, TIME, IX, = IX, F12.6 
18, F12.6 

First 
Solution 


Second 

Solution 


Last 
Solution 

where 

XXX = The TIME VALUE of the transient 

NODE = Stress point name 

TEMP = Temperature at this point 

This output is always produced. It can be reformatted with simple 
conversion programs to produce temperature input for almost any stress 
analysis code. 

An example of this output is in Appendix B. 


TIME = XXX 
NODE TEMP 


TIME = XXX 
NODE TEMP 
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5.2 NASTRAN TEMPERATURE FILE 


This optional output file is produced if the stress analysis program 
used is NASTRAN. The stress point temperatures are formatted so that they 
can be directly input into NASTRAN. The form of the file is: 


$ TIME = XXX 
TEMP NC NODE Temp 


Formats 

$, IX, TIME, IX, = IX, F12.6 
TEMP, 4X, 218, F8.2 


First Solution 


where 


Last Solution 


XXX = Time value of the solution 

NC = Case number (1 through number of solutions) 

NODE = Stress point name 

TEMP = Temperature of the stress point 

An example of this file is m Appendix B. 


5.3 MARC TEMPERATURE FILE 


This output file is produced if the stress analysis code to be used is 
MARC. The temperature information is generated as a set of Change State 
cards with a header card to indicate the time step from the heat transfer 
analysis. This data must be manually edited into the MARC input. 

Since MARC is a somewhat unique code m that it requires either centroid 
or integration point temperatures, an element connectivity file must also be 
included as input to TRANCITS. Currently, temperature input is available 
only for eight-noded stress solid elements (Element Type 7). 


5.4 HARD-COPY OUTPUT DESCRIPTION 


The hard-copy output of TRANCITS has basic information about the heat 
transfer model and the stress model to be interfaced. It starts with a 
TRANCITS header page, then it prints counts of the number of nodes and 
elements in the heat transfer model. It also prints the total number of 
solution time steps and the number of solutions processed for the run. 
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Next, it prints a header line saying "The following stress points were 
not contained by any heat transfer element" and that "SNOD XO YO 20 NEARE XN 
YN ZN DIST IWARN . " 

If no numbers appear under this header, all stress points in the stress 
analysis model wera inside a heat transfer element. Any stress point outside 
of the heat treat elements will appear under the header. The definitions of 
each of the columns are: 


SNOD 
XO \ 
YO 
Z0 ) 
NEARE 
XN 
YN 
ZN 

DIST 

IWARN 


Stress node name outside all heat transfer elements 


Original coordinates of this stress point 


Heat transfer element name nearest to this stress point 


Surface coordinates on nearest heat transfer element for which 
the temperature of this stress point was computed. 


Distance between (XO YO Z0) and (XN YN ZN) 

0 - Surfacing procedure converged normally 

1 - Surfacing procedure had difficulty m converging. 

(Check temperatures at this stress point.) 


NOTE: The two main points to look for in this printout are the magnitude 

of the distance parameter (if it is large, then the stress point 
was far outside the heat transfer model) and the IWARN parameter. 


The last portion of the hard-copy output deals only with finite- 
difference heat transfer codes. The header reads: 


"Corner Node Comparison For Time = XXX" 

"TCRITICAL For This Case Is Set to 99°" 

"NODE # MAXOT TBAR **C0MM0N H/T ELMS **" 

This printout is a consequence of the heat transfer element corner node 
calculation. Since the finite difference codes yield temperatures at the 
element controids and the element face centers, TRANCITS uses these tempera- 
tures to compute the element corner temperatures for each element. These 
corner temperatures are compared and if the difference of the temperature at 
the same corner is greater than 99°, then that corner is flagged in this 
printout. The column definitions are: 
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Node No. 

Max PT 
TBAR 

Common Heat Transfer 
Elements 


Heat transfer corner node name 

Maximum delta temperature at that corner 

Average temperature at that corner 

Heat transfer elements bounding the corner 


If this type of printout occurs it is an indication of large thermal 
gradients in the heat transfer model. If the MAXDT value is large there is a 
potential problem with the adequacy of the mesh density of the heat transfer 
mesh. When this situation occurs it is recommended that the heat transfer 
analyst be consulted about the validity of the results in the high gradient 
area . 


The last message printed is 

"THE INTERFACING PROCEDURE HAS BEEN COMPLETED." 
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6.0 RUN COMMANDS 


The current version of TRANCITS (Version 1.0) is operational on the NASA- 
Lewis IBM 370 system. 

Upon release of this user's TRANCITS manual, a formal procedure will not 
have been set up to run in production mode at NASA-Lewis. This procedure 
will be established in the near future and will be communicated to all users. 


21 



This Page Intentionally Left Blank 



APPENDIX A 


CONTENTS OF MARC POSTTAPE - VERSION K.l 


This appendix contains the documentation of the MARC Heat Transfer 
Input. Refer to the MARC manuals for more information. 
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APPENDIX A 


CONTENTS OF MARC POSTTAPE - VERSION K.l 


The following describes the contents of the 10 blocks found on the 
posttape. This information is the same for either the binary or the 
formatted posttape. 

BLOCK NUMBER DESCRIPTION 

BLOCK I TITLE 

NR;80A1 ; (TITLE (J) , J=1 , 70) TITLE (J) is Jth character 


BLOCK II 
NR;6I3;INUM 
LNUM 
MNUM 
NDEG 
NSTRES 


INOD 

IPSTCC 

IPSTYP 

NCRD 

NNODMX 


CONTROL INFORMATION 

Number of variables each element 

Number of nodes m mesh 

Number of elements in mesh 

Number of degrees of freedom each node 

Maximum number of integration points 
each element. 

NSTRES = 1 if ALL POINTS is not flagged 

NOTE: For constant strain elements 

having only one integration point 
NSTRES = 2 if ALL POINTS is flagged. 

Number of nodal variables - see note 
below. 

Connectivity and coordinate flag 
(0-none; 1-given) 

Type of tape (0 - binary, 1 - formatted 
2 - both) 

Number of coordinates each node 
Maximum number of nodes each element 
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IANTYP Analysis Type Flag: 

1 = Displacement - without reaction 

forces 

2 = Displacement - with reaction forces 

3 = Dynamics - without reaction forces 

4 = Dynamics - with reaction forces 

5 = Heat transfer 

6 = Rigid plastic flow - without reaction 

forces 

7 = Rigid plastic flow - with reaction 

forces 

8 = Harmonic analysis - without reaction 

forces 

9 = Harmonic analysis - with reaction 

forces 

10 = Harmonic analysis with dynamics - 

without reaction forces 

11 = Harmonic analysis with dynamics - 

with reaction forces 

12 = Heat transfer with fluxes 

13 = Joule heat 

14 = Bearing 

15 = Eigenvector (modal) 

16 = Eigenvector (buckle) 


ICOMPL 

Set to 0 if real analysis 
Set to 1 if complex analysis 

NBCTRA 

Number of nodes with transformations 

IPOSTR 

Posttape revision number, 1 for this 
release (J3) 

IDM4 



I DM5 

I 

Not used; reserved for future 

IDM6 

IDM7 

j 

expansion 
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Note : 


INOD = {NDEG }* JNODE 

IF IANTYP = 1; 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 


JNODE = 1 
on 2 

3 

4 
1 
1 
2 
1 
2 

3 

4 
2 
3 

5 
1 
1 


IANTYP 
IANTYP 
If IANTYP 
If IANTYP 


= 15 only appears during 
= 16 only appears during 
= 8 and ICOMPL = 1 JNODE 
= 9 and ICOMPL = 1 JNODE 


subincrements 
subincrements 
= 2 
= 4 


BLOCK III 


NR;6I3 ; JPLOT(J) 


CODE NUMBERS ASSOCIATED WITH ELEMENT 
VARIABLES 

INUM records; J=1 , INUM 

Element variable code (see MARC Manual 
Volume C, p C3.10-7) 


BLOCK IV 


NR;6I3;LM(1) 

LM(2) 

LM(3) 

LM(NN0DMX+2) 


CONNECTIVITY LIST 

If IPSTCC is zero this block will be 
omitted 

MN1JM records; J=1 ,MNUM 
Element type 

Number of nodes in this element 
First node of Element J 
Last node of Element J 


BLOCK V 


COORDINATE LIST 


If IPSTCC is zero this block will be 
omitted . 
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NR;6E13.6; SUM(l) 
SUM(2) 
SUM(NCRD) 


BLOCK VI 


NR; 61 13; LM(I) 


BLOCK VII 


NR;6E13.6; D(l,l) 
D(2,l) 
D(3 , 1) 
D(l,2) 
D(2,2) 
D(3,2) 
D(l,3) 
D(2,3) 
D(3,3) 


Blocks VIII, IX, and 
BLOCK VIII 


NR;6E13.6 Xl(l) 
Xl(2) 


XI (3) 


LNUM Records; J=1,LNUM 
First coordinate of Jth node 


Last coordinate of Jth node 


TRANSFORMATION LIST 


If NBCTRA is zero or if IPSTCC is zero, 
this block will be omitted. 

Binary 1 record 

Formatted (NBCTRA-l)/6 + 1 records 

List of nodes which have transformations 
applied . 


TRANSFORMATION - DIRECTION COSINES 

If NBCTRA is zero this block will be 
omitted. 

One record for each node listed in 
Block VI if binary tape. 

Two records each node if formatted tape. 
Transformations are for local to global 
system. 


are repeated for each increment. 

INCREMENT, TIME, AND FREQUENCY 
Transient Time 

Increment number is a real formed 
as I + J/100 

I is the static increment number 
J is either the harmonic subincrement 
or the Eigenvector number. 

Frequency 
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XI (4) 


Flag to read new Blocks II, III, IV, V, 
VI, VII. Set to 1 to read these blocks 
again. This is only used in a rezoning 
analysis . 

XI (5) IANTYP analysis type flag 

XI (6) Not used; reserved for future expansion. 


BLOCK IX VALUES OF ELEMENT VARIABLES 

If INUM is zero, this block is omitted. 

If IANTYP=15 or IANTYP=16, this block 
is omitted. 

MNUM*NSTRES Records 

DO 411 11=1 ,MNUM 
DO 413 JJ=1 ,NSTRES 

NR;6E13.6; SUM(l) Value of 1st element variable for 

• Element II, Integration Point JJ. 


SUM(INUM) Value of INUM^ element variable for 

Element II, Integration Point JJ. 

413 CONTINUE 
411 CONTINUE 


BLOCK X 


NR;6E13.6; SUM(l) 

SUM(NDEG) 

SUM(NDEG+1) 

SUM(2*NDEG) 

SUM(INOD) 


VALUE OF NODAL VARIABLES 

If INOD is zero, this block is omitted. 
LNUM records 

Nodal displacements, velocities, 
accelerations, and reactions 


During static or dynamic increments : 


The first NDEG quantities are: 


IANTYP = 1, 2, 3, 4, 8, 9, 10, 11 
Total nodal displacements 
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IANTYP = 6, 7 

Total nodal velocity 

IANTYP =5, 12, 13, 

Total nodal temperature 

IANTYP = 14 

Pressure 

The second NDEG quantities are: 

IANTYP =2, 7, 9 

Nodal reactions at boundary 
conditions; residual load 
correction elsewhere 

IANTYP = 3, 4, 10, 11 

Nodal velocity 

IANTYP =12, 13 
Fluxes 

IANTYP = 14 

FORCEX 

The third NDEG quantities are: 

IANTYP = 3, 4, 10, 11 

Nodal accelerations 

IANTYP =13 

Nodal voltages 

IANTYP =14 

FORCEY 

The fourth NDEG quantities are: 

IANTYP =4, 11 

Nodal reactions at boundary 
conditions; residual load 
correction elsewhere 

IANTYP =14 

FORCE Z 

The fifth NDEG quantities are: 

IANTYP =14 

MASS FLUX 
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During subincrements : 

The first NDEG quantities are: 

IANTYP = 15 

Nodal components of dynamic 
mode 

IANTYP =16 

Nodal components of buckle 
Eigenvector 

IANTYP =8, 9, 10, 11 

Real components of harmonic 
displacement 

The second NDEG quantities are: 

IANTYP =9, 11, and ICOMPL = 0 

Real components of harmonic 
reaction force 

IANTYP =8, 9, 10, 11, and ICOMPL = 1 
Imaginary components of 
harmonic displacement 

The third NDEG quantities are: 

IANTYP =9, 11, and ICOMPL = 1 

Real components of harmonic 
reaction force 

The fourth NDEG quantities are: 

IANTYP =9, 11, and ICOMPL = 1 

Imaginary components of 
reaction force 

NOTE: There will always be INOD quan- 

tities; the record may be padded 
with zeroes m subincrements. 

IF IANTYP = 14, same as INCREMENTS. 
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NOTES: 


NR Indicates the beginning of a new record 

for the binary posttape. 

NR;Format Indicates the beginning of a new set of 
information which is to be read with the 
following format. 
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APPENDIX B 


This Appendix shows examples of the following file formats: 

• Neutral Heat Transfer Input 

• SINDA Temperature Input 

• SINDA Geometry Input 

• MARC Postinput 

• Stress Model Coordinate Input 

• Stress Model Element Input 

• User Control Input 

• Neutral Temperature Output 

• NASTRAN Temperature Output 
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NEUTRAL HEAT TRANSFER INPUT FILE 


This is the Test Case data. This is the new file form to be 
delivered for NASA's version — we are using a very simple test 
case, as follows: 


09 - 15-83 (SR MARTIN ) 


0 5 ^ 20 1 20 11 1 DIRECTORY 


1 

1.000000 

1 . 000000 

1 . 000000 




i - 

1 . 000000 

2.000000 

1.000000 




3 

1 .000000 

3.000000 

1.000000 




11 

2.000000 

1 .000000 

1.000000 





> 

w 




NODAL 

243 

5.000000 

3 . 000000 

3.000000 




251 

6 . 000000 

1 . 000000 

3 . 000000 




25Z 

6 . 000000 

2.000000 

3.000000 




253 

6 . 000000 

3.000000 

3.000000 




1 

1 2 

12 1 1 

101 102 

112 

111 


2 

2 3 

13 12 

102 103 

113 

112 


11 

11 12 

22 21 

111 112 

122 

121 

ELEMENT 

12 

12 13 

23 22 

112 113 

123 

122 

CONDUCTIVITY 

21 

21 22 

32 31 

121 122 

132 

131 


22 

22 23 

33 32 

122 123 

133 

132 


31 

31 32 

42 41 

131 132 

142 

141 



* 


141 

141 142 

152 

151 

241 

242 

252 

251 


142 

1 

L. 

142 143 

123456 
123456 

153 

152 

242 

243 

253 

252 

ELEMENT 

FACE 


11 123456 

12 123456 

21 123456 

22 123456 


t 


132 

123 

456 


141 

123456 


142 

123 

0 . 

456 


1 

200.0 

100.0 

200.0 

2 

200.0 

100.0 

200.0 

11 

400 . 0 

300.0 

400.0 

12 

400.0 

300.0 

400.0 

21 

600.0 

500.0 

600.0 

22 

600.0 

500.0 

600.0 

31 

800.0 

700 . 0 

300 . 0 


300.0 

200.0 

200.0 

200.0 

300.0 

200.0 

200.0 

200.0 

500.0 

400.0 

400.0 

400.0 

500.0 

£ 00.0 

400 . 0 

400.0 

700.0 

600 . 0 

600.0 

600.0 

700 . 0 

600.0 

600.0 

600 . 0 

900.0 

800.0 

800 . 0 

900 . 0 


TEMPERATURE 
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SINDA TEMPERATURE FILE 


255 


101 

102 

103 

1?4 

195 

136 

107 

138 

139 

110 

11! 

112 

113 

111 

115 

116 

117 

118 

119 

123 

121 

122 

123 

124 

125 

126 

127 

128 

129 

133 

131 

IX 

133 

134 

135 

136 

137 

138 

139 

145 

141 

142 

143 

144 

145 

146 

147 

148 

149 

155 

151 

152 

153 

154 

155 

156 

157 

153 

159 

165 

161 

162 

163 

164 

165 

166 

167 

168 

169 

173 

171 

172 

173 

174 

175 

176 

177 

178 

179 

153 

181 

182 

183 

184 

185 

291 

252 

233 

234 

235 

296 

237 

233 

2 89 

219 

211 

212 

213 

214 

215 

216 

217 

218 

219 

229 

221 

222 

223 

224 

225 

226 

227 

223 

229 

229 

231 

232 

233 

234 

235 

236 

237 

238 

239 

240 

241 

242 

243 

244 

245 

246 

247 

248 

249 

250 

251 

252 

253 

254 

255 

256 

257 

258 

259 

260 

261 

262 

263 

264 

265 



371 

372 

373 

374 

375 

376 

377 

378 

379 

380 

381 

382 

283 

384 

385 







0.39165E 0* 

0.59459E 03 5.49561E 03 0.53582E 33 0.49601E 03 0.50756E 53 0.49754E 03 

0.50943E 92 S.49779E 83 0.511*71 03 3.5001K 93 9.59953E 03 S.49697E 03 

0.50257E 83 0.48327E 03 3.50192E 03 0.49241E 33 0.48952E 83 9.49366E 83 

0.48742E 03 0.49I50E 03 0.*8658E 93 I.48B96E 03 S.48728E 93 0.4ES89E 03 

0.48962E 03 0.48934E 13 I.49932E 83 0.48958E 03 0.49924E 83 0.48975E J3 

I.48925E 03 9.48939E 13 0.48886E 93 8.4393« 03 0.48949E 03 S.*8881E 03 

I.49317E 03 0.49163E 03 9.49448E 03 0.49218E 83 0.49580E 83 3.49246E 83 

8.49543E 03 I.49286E 83 M9631E 83 S.49391E 83 I.47638E 83 J.47645E 03 

I. 47225E 03 8.49114E 93 0.48905c 03 0.4847*£ 13 1.49605c 93 8.49548E 93 

9.49394E 03 0.49825E 13 0.49749E 83 I.49597E 53 M9986E 03 0.4986S 03 

J. 49692E S3 0.50156E 13 0.40935c 03 0.49653E 03 0.50273E 03 0.49954E 03 

0.49590E 03 0.50356E 03 0.49991E 13 0.49590E 03 0.50445c 03 M9964E 03 


LIST OF 
ELEMENTS 


LIST OF 
TEMPERATURES 



SINDA GEOMETRY FILE INPUT 


SINDA GEOMETRY 

FILE - 3D 

COMS'JSTOR LINER 

MODEL 

500 255 

8 



1 

0 23640 

15 75804 

0 . 

2 


1 5 792 JO 

0 . 

3 

0 53550 

15 74303 

0 . 

4 

0 53730 

15,77885 

0 . 

5 

0 29100 

15 82606 

0 . 

6 

0 _ o401 0 

16 61478 

0 . 

7 

0.78450 

15 72802 

0 . 

8 

0.78690 

15 76570 

0 . 

9 

0.78930 

15 80339 

0 . 

10 


J~"7 

0 . 

1 1 

1 . 03600 

15.75256 

0 . 

12 

1 . 03840 

15.79210 

0 . 

13 

1 . 10760 

15 70856 

0 

14 



0 . 

15 

1.11410 

15 81757 

0 

16 

1 . 16370 

15 70524 

0 

17 

1 16690 

15.75970 

0 . 


t t t t 


717 



R" 






718 


2.11210 

14 

93652 


5 43645 



71 9 


2. 11230 

14 

91201 


5 42753 



720 


2 11250 

14 

68749 


5.41860 



721 


■ 1 ■■ 1 II — 







722 


2. 19960 

14 

93681 


5.43656 



723 


2 19980 

14 

91239 


5.42767 



724 


2 20000 

14 

88807 


5 41882 



725 


2 20020 

14 

663q5 


5 40993 



101 

203 

204 

202 

201 

3 

4 

2 

1 

1 23456 

102 

204 

206 

205 

202 

4 

6 

5 

2 

123456 

103 

207 

208 

204 

203 

7 

8 

4 

3 

123456 

1 04 

208 



533" 

8 

9 

6 

4 

1 23456 

105 

210 

21 1 

208 

207 

10 

1 1 

8 

7 

1 23456 

106 

21 1 

212 

203 

208 

1 1 

1 2 

9 

8 

1 23456 

107 

213 

214 

21 1 

210 

13 

14 

1 1 

10 

12345S 

108 

214 

215 

"5T2 

mi 


■a 

gmp 

1 1 

1 23455 

109 

216 

217 

214 



Kfl 


13 

123456 

1 10 

217 

21 8 

215 



S£I 


14 

1 23456 

1 1 1 

221 

222 

217 



M'J'M 


16 

123456 

1 1 2 

222 


218 

217 

55 

5 IT 

1 8 

1 7 

123456 

1 13 

228 

229 

222 

221 

28 

29 

22 

21 

1 23456 

1 14 

229 

230 

223 

222 

29 

30 

23 

22 

123456 

115 

234 

235 

229 

229 

34 

35 

29 

28 

1 23456 

1 16 

235 

235 


m'4'41 

35 


30 

29 

123456 

117 

224 

220 

218 

223 

24 

20 

1 8 

23 

1 23456 

118 

220 

219 

215 

218 

20 

19 

15 

18 

123456 

119 

227 

226 

220 

224 

27 

26 

20 

24 

1 23456 

120 

226 


219 

220 

26 

25 

19 

20 

123456 

121 

233 

232 

226 

227 

33 

32 

26 

27 

123456 

122 

232 

231 

225 

226 

32 

31 

25 

26 

123456 

123 

239 

238 

232 

233 

39 

38 

32 

33 

123456 

124 

238 

537 

231 

232 

38 

37 

31 

32 

123856 

125 

238 

239 

242 

241 

38 

39 

42 

41 

123455 

126 

237 

238 

241 

240 

37 

38 

41 

40 

123456 

127 

241 

242 

245 

24 4 

41 

42 

45 

44 

1 23456 

128 

240 

241 

244 

243 

40 

41 

44 

43 

1 23856 



^ ^ ^ ^ 
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MARC POSTFILE/TAPE INPUT 


C0M3USTSR MARC POST FILE SIMULAT1CS 


0=00255 1 

ET 

T 

113 6 

0 1 

5 

0 

43 8 3 4 

2 

1 

1261291 2TT?6 
43 8 4 6 

129131 1 33 1 27 

5 

2 

43 8 7 8 

4 

3 

1 321 33 1 291 2o 
43 8 8 9 

133134131 129 

6 

4 

43 81011 

8 

7 

‘ 135136133132 

43 8 11 12 

136137134133 

9 

6 

43 3 13 14 

1 1 

10 


* 


495494490491 
43 6369368354365 

494493489430 

43 63753743' 7 0371 

"500499495435 
43 8374373369370 

499493494495 

43 6373372368369 

496497493494 

0 286400E 00 0 1 57580E 02 0 

0 28S700E 00 0 1 57920E 02 0 

0 535500E 00 0 1 57430E 02 0 

0 537800E 00 0 1 57789E 02 0. 

0.291 OCOE 00 0 1 58261 E 02 0. 
0 540100E OC 0 158148E 02 0. 

0 764500E 00 0 157260E 02 0 

* 


0 


~3T" 

tr 


ur 

TT 

"542753E' 

01 

0 

21 1250E 

01 

0 

146875E 

02 

0 

541860E 

01 

0 

21 1270E 

01 

0 

148631 E 

02 

0 

540972E 

01 

0 

219960E 

01 

0 

149368E 

02 

0 

543656E 

01 

0 


TTT 

ir 

1 49 ’ 24t 

K!!H 

~c 

5427&7E 

01 

0 

2200C0E 

01 

0 

1 46881 E 

02 

0 

541682E 

01 

0 

220020E 

01 

0 

148637E 

02 

0 

540393E 

01 

0 

393700E 

04 

0 

1000CCE 

01 

0 


0 0 300000E O’ 










0 

499432E 

03 







0 

512448E 

03 







0 

500855E 

03 







0 









0 

491059E 

03 







0 

509778E 

03 







0 

. ... . , ' * ' 

49651 8E 

03 

— 








0.510112E 03 
0.494647E 03 
0.482880E 03 

* 
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STRESS MODEL COORDINATE INPUT FILE 


1 

1 34980 

15 86530 

0. 

2 

1 . 34960 

li; _ a3SoO 

U. 

3 

1 .34950 

15 60435 

1 . 03S43 

4 

1 34930 

15 83179 

1 . 03823 

5 

1 .39740 

15 80455 

1 03644 

5 

T73S730 


1 03822 

7 

1.39740 

15 83650 

0. 

0 

1 .39753 

15 66570 

0 . 

9 

1 . 34920 

15 81080 

0. 

TO 


T5 77F91 

1 . 03463 

1 1 

1 39727 

15.77731 

1 . 03465 

12 

1 .39727 

15 81120 

0. 

1 3 

1 . 34920 

15 67554 

2 06372 

14 




15 

1 34960 

15 73C07 

2 07090 

16 

1 39727 

15 67593 

2 06378 

17 

1 39740 

15 70300 

2 C6734 

TS 


15' 72997 

2 07069 

1 9 

1 44533 

15 67633 

2 06363 

20 

1 44533 

15 77771 

1 03463 

21 

1 44533 

15 61160 

0 

22 



2 06735 

23 

1 .44527 

15.72997 

2 07089 

24 

1 .44530 

15 80465 

1 .03645 

25 

1 .44527 

15 83169 

1 . 03822 

26 



0 . 

27 

1 .44527 

15 86570 

0 

28 

1 .34920 

15 5071 1 

3.03399 

29 

1 34950 

15 53408 

3 08935 

30 

■Ec££X-I4IEI 


3 09472 

31 

1 .39727 

15 50750 

3.08407 

32 

1 39740 

15 53428 

3.08940 

33 

1 39753 

15 56095 

3 09470 

34 


■usucai 

3 08415 

35 

1 .44530 

15 53437 

3 08942 

36 

1 .44527 

15.56C95 

3 09470 

37 

1 .49340 

15 50628 

3 03423 

38 



3 08945 

39 

1 .49300 

15 56C86 

3 09463 

40 

1 49340 

15 67673 

2 06386 

41 

1 49340 

15 7761 1 

1 03471 

42 


— 15 61 200 

0 . 

43 

1 .49320 

15 70330 

2.06738 

44 

1 49300 

15 72987 

2 07088 

45 

1 49320 

15 80435 

1 03646 

45 



0 . 

47 

1 .49300 

15 83159 

1 .03821 

48 

1 .49300 

15 86560 

0 . 

49 

1 34920 

15 27206 

4 09214 

50 



4 09925 

51 

1 .34980 

15 32519 

4. 10637 

52 

1 .39727 

15 27245 

4.09224 

53 

1 . 39740 

15 29682 

4 09931 
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STRESS MODEL ELEMENT INPUT FILE 


1 

1 1 452 


1 

790 

789 

7S6 

788 

787 

785 

783 

784 

2 

7S9 

782 

781 

78G 

785 

780 

779 

783 

3 

762 



/61 


/b i 

/bU 

779 

4 

763 

726 

725 

762 

751 

724 

723 

760 

5 

76 6 

786 

777 

776 

784 

783 

775 

776 

6 

786 

7ei 

770 

111 

783 

779 

769 

775 

7 

761 

762 

'~759 

BffcB 

/ /9 

/bU 

/5o 

/t>9 

6 

762 

725 

722 

759 

760 

723 

721 

758 

9 

778 

777 

753 

754 

776 

775 

751 

752 

10 

111 

770 

750 

753 

775 

769 

749 

751 

1 1 


/oy 

/Jb 





775 

12 

759 

722 

717 

736 

758 

721 

716 

735 

1 3 

78 -7 

785 

783 

784 

774 

773 

771 

772 

14 

785 

780 

779 

783 

773 

766 

767 

771 

1 5 


751 

/oQ 

/ /9 

/bO 

7b / 

756" 

757 

1 6 

761 

724 

723 

760 

757 

720 

719 

756 

17 

784 

783 

775 

lie 

772 

11' 

765 

766 

1 6 

733 

779 

769 

775 

771 

151 

764 

765 

1 9 

779 

7bU 



/ = / 

/5t> 

/So 

757 

20 

760 

723 

721 

758 

756 

719 

718 

755 

21 

776 

775 

751 

752 

766 

765 

747 

74 8 

22 

775 

769 

749 

751 

765 

764 

746 

747 


t 


427 

497 

495 

492 

^53 

457 ~ 

456 

452 

^55 

428 

496 

495 

491 

492 

456 

455 

451 

452 

429 

658 

638 

637 

657 

630 

629 

627 

628 

430 

638 

60S 

607 

637 

629 

60C 

599 

627 

431 


07 i 




5bC 

s=y 

" 599 ‘ 

432 

657 

637 

636 

656 

628 

627 

625 

626 

433 

637 

607 

60S 

636 

627 

599 

596 

625 

434 

607 

570 

569 

606 

699 

559 

558 

596 


■$f4i 

63b 

635 

65= 

626 

625 

623 

S"24 

436 

636 

606 

605 

635 

625 

596 

597 

623 

437 

606 

569 

565 

605 

598 

556 

557 

597 

438 

655 

635 

631 

651 

624 

623 

621 

622 

439 

635 

60S 

601 

651 

623 

597 

596 

621 

440 

6C5 

565 

561 

601 

597 

557 

556 

596 

441 

630 

629 

627 

628 

595 

594 

592 

593 

442 

629 

600 

599 

627 

594 

591 

590 

592 

443 




555 

551 

555 

5=4 

590 

444 

628 

627 

625 

626 

593 

592 

588 

569 

445 

627 

599 

598 

625 

592 

590 

587 

538 

446 

599 

559 

558 

598 

590 

554 

553 

587 

447“' 


■444 


■4^1 


56b 

58 = 

586 

448 

625 

598 

597 

623 

568 

587 

584 

585 

449 

598 

558 

557 

597 

587 

553 

552 

584 

450 

624 

623 

621 

622 

586 

585 

582 

583 

451 

623 

597 

"596 

621 



56 1 


452 

597 

557 

556 

596 

584 

552 

551 

581 
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USER CONTROL FILE 


SZMAF 

ISPOC=1 , 

NTMD0=3, 

NUM2AP/NZA°= 1 ,-11,-2,51, 
nTST r/ ro ' , ' 30 ' , si , ' 

* 
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NEUTRAL TEMPERATURE OUTPUT FILE 


TIME = 3916 030000 


1 

775 ~ 

251T75 

2 

475 

362076 

3 

473 

0S1 939 

4 

472 

954633 

5 

477 

105648 

6 

477 

014626 

7 

479 

731720 

e 

476 

416409 

9 

776" 

380 043 ~ 

10 

476 

101242 

1 1 

460 

371734 

12 

481 

297924 

1 3 

"477 

1 93 

14 

475 

659309 

15 

473 

984970 

16 

482 

275578 

1 7 

479 

5752S6 

18 

477 

891411 

19 

487 

51 2489 

20 

485 

714573 

TIME s 

3947 000000 

1 

913 

990150 

2 

914 

288366 

3 

912 

6735S2 

A 

5T2" 

"V20973 — 

5 

921 

552902 

6 

921 

516510 

7 

925 

217186 

8 

923" 

1)2 1 370 

9 

914 

72071 1 

10 

912 

963562 

1 1 

923 

802025 

1 2 

926 

U4 39 » 5 

13 

91 1 

039658 

14 

910 

560577 

15 

910 

516815 

T 5 

922 

a i 4a63 

17 

91 9 

726219 

18 

918 

977898 

19 

937 

8C2277 

20 

33/ 

246/73 

TIME = 

4055 000000 

1 

1027 

800430 

2 

1028 

4503 1 7 

3 

1 028 

727051 

4 

1023 

784714 

5 

1044 

29901 1 

6 

1044 

251038 

7 

1 045 

390854 

8 

T 043 

"329732 

9 

1027 

9991 15 

10 

1029 

331284 

1 1 

1045 

346954 

12 

l 0^6 

466S4 1 

13 

1030.756336 

14 

1030.794373 

15 

1030 

31 96e7 

1 6 

1 U4 9 

595552 

17 

1046 

791626 

1 8 

1 045 

744339 

19 

1077 

498398 

20 

1D74' 

42 1 733 


SOLUTION 1 


SOLUTION 2 


SOLUTION 3 
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NAS TRAN TEMPERATURE FILE 


S TIME = 391 o OCCOOO 



1 T 

— 4 72" 2 b 

TEMP 

1 2 

475 36 

TEMP 

1 3 

473 06 

TEM P 

1 4 

472 96 

TEMP 

1 5 

477' T1 

TEMP 

1 6 

477.01 

TEMP 

1 7 

479 73 

TEMP 

i a 

476 42 

TER? 

1 3 

4 75 36 

TEMP 

1 10 

476 10 

TEMP 

1 1 1 

480 37 

TEMP 

1 12 

481 30 


1 1 3 

'477 67 

TEMP 

1 14 

475 66 

TEMP 

1 15 

473 98 

TEMP 

1 15 

482 26 

TEMP 

1 17 

'479 58 

TEMP 

1 18 

477 89 

TEMP 

1 19 

487.51 

TEMP 

1 20 

485 71 

TEHP 

1 2 1 

437 03 

TEMP 

1 22 

484 75 

TEMP 

1 23 

483 24 

TEMP 

1 24 

482 65 

TEMP 

1 25 


TEMP 

1 26 

485 24 

TEMP 

1 27 

462 10 

TEMP 

1 28 

479 20 


I ZS~ 

4 76.43 

TEMP 

1 30 

476 07 

TEMP 

1 31 

483 04 

TEMP 

1 32 

480 08 


1 

479 67 

TEMP 

1 34 

487 76 

TEMP 

1 35 

484 . 96 

TEMP 

1 36 

484 60 
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APPENDIX C 


SINDA FORTRAN OUTPUT SUBROUTINES 


This appendix contains listings of the SINDA Fortran output subroutines 
used to produce the SINDA Temperature Input File. Refer to the SINDA Manual 
for proper use of these routines. 
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SINDA FORTRAN OUTPUT SUBROUTINE 


C 

c 


BCD SEmfcut 1 oh 1 


ASSIGN DYNAMIC STORAGE PECUIREMEHTS 
SIZE OF X DEPENDS ON NUMBER OF NODES, SOLUTION ROUTINE. 
AND SUBRGU7 1NES CALLED 

• 

COMMON/L QGI C/LNODE . LCCND , L CONST , LARK AY 
DIMENSION XC200 ) 

ND1M=20D 
NTH =0 


THE FOLLOWING CODE FILLS THE DYNAMIC STORAGE AP.PAY *X’ (‘NX*) 

WITH THE NODE NUMBERS IN THE SAME CRDSP AS THE »T» ARP AY . 

TriEH BEGINNING WITH THE FIRST NODE HUMBER. THE NODES AP.E , 

WRITTEN TO THE FILE USED BY ’TEMPSCAN’ #r J. IM _ ^ / . 

S5S REFER TO SINDA USER’S MANUAL. PAPF S-S «$$ v 1 ‘ X 

icji'S Jr*’/ Ue meu'£,‘ 

NOTE: THIS HAPPENS ONLY ON THE FIRST CALL TO EXECUTION DUE 

TO THE CHECK ON THE STATUS OF ’LNODE' . 

IFCLNODE.GE.O) GO TO ZOC f? 

IOC = NTH '* 

LEN = NDIM 
»LL HNREADf! 


SJC . zeniths, 
-c>n ly-..c*lled y one*?. 


WRIT EC 23,100)NNT - 

ICO FONUTCU) ~ r L%r° l " r 

WRITEC23.I10) (NXtI*NDIM),I = l.NNT) ^rehe* 'fa i>U 

110 FORMAT (1018) ^ 

ioo continue feph+e* omv 


c 

c 

c 

c 

c 

c 

F 

F 

F 

F 

C 

c 

C 

C 

C 

c 

c 

c 

c 

F 

F 

F 

F 

C 

F 

C 

F 

c 

F 

F 

C 

C 

F 


EXECUTE PPCGPAM 
CNFFUID 

cscnrii’ 

END 


C 


c 

r 


BCD 3VAPIABLES 2 
END 


BCD B OntPUT 


TPPNTF 


W° 1 7 E NODE TEMPERATURES TO TAPE 
KTE5 T =KTESTM 
W'PITE(6 . 130 )KTEST. TIMES 

130 FOFMAT t SX , "S S $ NOTE- DATA SET =",15." WRITTEN TO TEMPSCAN", 

* " FILE AT TIME = " , FI 0 3) 

WPITE ( 2 3 . 1 A 0 ) TIMES WtTTf ~ r v & If 


1 A 0 


i:?ITE(23.1A0)CTCI),I = l,NST) v 
FORMAT ( BE) 2 . S ) 

END 


>t* yre> ot > v* 


, r 

’r+mmr* 


-fr 


'flip. 


c 

c 

c 

c 

c 


c 

c 

F 

F 

F 

F 

F 

F 

F 

C 
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C 

c 



End of Document 



